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\0 1 Abstract 

o 

It is known that some of the deeply-virtual-Compton-scattering observables, for instance, 
the beam-charge asymmetry in the hard electroproduction of real photons on the nucleon, 
are extremely sensitive to the magnitude of D-term appearing in the parameterization of the 



Oh! 



generalized parton distributions. We report a theoretical analysis of both the isoscalar and 
isovector parts of the nucleon D-term within the framework of the chiral quark soliton model, 
without recourse to the derivative expansion type approximation used in previous works. 
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^ ■ Growing attention has recently been paid to the studies of the so-called generalized parton 

distributions (GPDs) not only because they provide us with the unique means to experimen- 
tally access the quark orbital angular momentum in the nucleon but also because they offer 
the most detailed information on the underlying quark-gluon structure of the nucleon [1] - [7] . 
However, since the GPDs are functions of three kinematical variables and since they appear as 
complicated convolution integrals in the cross section formulas of the deeply-virtual Compton 
scatterings (DVCS), deeply- virtual meson productions (DVMP) etc., a suitable parameteri- 
zation of them is practically unavoidable. The most popular parameterization of GPDs is to 
use the double distributions p],[l] supplemented with the so-called D-term [8]. It turned out 
that some of the DVCS observables, for instance, the beam-charge asymmetry in the hard 
electroproductions of real photons on the nucleon, are extremely sensitive to the magnitude of 
the D-term necessary in the above parameterization [5], [9]. 

So far, there has been only a limited number of theoretical studies on the nucleon D-term. 
The first estimate of the D-term is based on the chiral quark soliton model (CQSM), more 
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precisely, on the CQSM predictions for the GPD H u+d (x,^,t) given by Petrov et al. under 
the derivative-expansion type approximation [TO]. Using these predictions, the authors of [5] 
as well as of jTT] estimated several Mellin moments of H u+d (x,^,t) at various values of £ and 
t and extrapolated them to t = 0. The results are then fitted to the Gegenbauer expansion 
given as 

oo 

D«(z) = (1-z 2 ) d nCl /2 (z). (1) 

q=u,d,s,--- n=l,odd 

This led them to the following estimate for the expansion coefficients : 

d x ~ -4.0, d 3 ~-1.2, d 3 ~ -0.4, (2) 

with higher coefficients being small. We recall here that the D q (z) in the flavor- singlet channel 
mixes under evolution with the corresponding gluon D-term D 9 (z). The numbers quoted 
in Eq.@ corresponds to the values at a few GeV scale. At the model energy scale around 
600 MeV, their result corresponds to (see the footnote of [12]) 

d x ~ -8.0. (3) 

Although used in many recent phenomenological analyses of the DVCS and the DVMP pro- 
cesses, such an estimate of the D-term coefficients is of highly qualitative nature. A little more 
direct estimate of d±, the first coefficient of the Gegenbauer expansion, was made by Schweitzer 
et al. within the same model [12]. By starting with the model expression for the unpolarized 
GPD H(x,£,t), they derived a closed formula for d%. They also estimate its numerical value 
by using the derivative-expansion type approximation to find that 

d 1 ~ -9.46, (4) 

at the model energy scale. After a simple estimate of the scale dependence, the author of [12] 
got the number 

di ~ -4.7 at a few GeV 2 , (5) 

which they claim is qualitatively consistent with the results of [H],[5]. Although of preliminary 
nature, there also exists a lattice QCD study of the generalized form factors Cf (£) [T3], forward 
limit of which are related to d±. 

Clearly, in spite of its phenomenological importance, we must say that our knowledge on the 
precise magnitude of the D-term is still rather poor and uncertain. In view of the circumstance 
above, we think it useful to evaluate the most important parameter of the D-term, i.e. the first 
coefficient of the Gegenbauer expansion of the D-term within the framework of the CQSM, 
without recourse to the derivative-expansion type approximation. In the present paper, we try 
to estimate not only the isoscalar part of d\ but also its isovector part, which is sub leading in 
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the 1/N C expansion. For the sake of comparison, we also derive the theoretical expression of 
d\ +d in the familiar MIT bag model, and see what prediction it gives. 

To fix the normalization convention of the relevant quantities, it would be convenient to 
start our investigation with the nucleon matrix elements of the quark and gluon parts of the 
(symmetric) QCD energy-momentum tensor parameterized by four form factors as |14j.[2] 



(p'lWO) \p) 



+ 



N(p') 
d q ' 9 (t) 



5M 



N 



N(p), 



(6) 



where P = (p' + p)/2, A = p' — p, and t = A 2 . Here f q = qy^ V „}g is the QCD energy- 
momentum tensor of the quark with flavor q, while = G a ^ a G a av + \ G 2 is the corresponding 
gluon part. (The form factor c(t) accounts for nonconservation of the separate quark and gluon 
parts of the energy-momentum tensor. They must satisfy the constraint c q (t) + c g (t) = 
due to the conservation of the total (quark plus gluon) energy-momentum tensor.) The form 
factors in Eq.fjH]) are related to the 2nd Mellin moments of the familiar unpolarized GPDs 
H q (x,£,t) and E(x,£,t) as 



J xH q (x,£,t)dx 



'-1 
-i 



M?{t) + -d9{t)e, 
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q 



J xE q {x,£,t)dx = 2J Q {t) - M 2 Q (t) 



V(tK 2 . 



(7) 
(8) 



Here, the suffix Q denotes the summation over all quark flavors, for example, J®(t) = 
J2 q =uds - - J q (t)- (Practically, we confine here to the light-quark components of two flavors, 
which means that Q = u+d.) The sum of the above two equations with t — gives the famous 
Ji's angular momentum sum rule [2],[3] : 



^2 j x{H q (x,{,0) + E q (x,£,0))dx = 2J Q , 



(9) 



with J® being the total angular momentum carried by the quark fields in the nucleon. 

The interest of our present study is the forward limit of d®(t), which just corresponds to 
the first coefficients in the Gegenbauer expansion of the so-called D-term, i.e. 



d\ 



u+d 



d Q {0). 



According to Polyakov [TJ], the constants d^ +d can be expressed as 

d 3 rT$(r 



ru+d _ M n 

dl ~ 2 



3 



(10) 



Here T® u (r) is the net quark contribution to the static energy- momentum tensor density of 
the nucleon defined by 



7?» = ^ / ^e- lA - r (p,S\f^(0)\p,S), (12) 



J_ f d*A lA . r 
^ v ' ; 2M W J (2vr) 3 ^>"i-^ 

with T® u (0) the quark part of the QCD energy-momentum tensor operator. As expected, the 
D-term is seen to contain valuable information on the distribution of energy-momentum tensor 
inside the nucleon. 

As a warm-up, let us first evaluate d\ in a simple model of baryons, i.e. the MIT bag 
model. We must calculate 

A/f 1 
< d (MIT) = --^{m gs \r 2 [oc.r^ -r - -cx-^}\S> gs ), (13) 

with a % = 7°7* the standard Dirac matrices. Here 

jJkr) I (I = 0) j = 1/2, m) 
*M = N\ Jo{ 1 1 1 . ' \ , (14) 
ijiikr) | (/ = 1) j = 1/2, m) 

with 

iV~ 2 = 2R 2 (kR-l)j (kR), uj = kR ~ 2.043 (15) 

is the ground state wave function with i? being the bag radius. After some manipulation, we 
easily find that 

< +d (MIT) = -M N N 2 k | jf j^kr) r A j^kr) dr - j (kr) r 4 j 2 (jfcr) dr J . (16) 

For an order of magnitude estimate for <i" +<i (MIT), here we use the bag model parameter 
adopted by Jaffe and Ji [T5], i.e. M N R ~ 4.0 c^o- Eq. lFIBl) then gives 

^ +d (MIT) ~ -0.716, (17) 

which turns out to be many times smaller in magnitude than the previous estimates based on 
the CQSM, d u + d ~ - (8.0 ~ 9.5). 

Next, we derive the theoretical expression for d\ within the CQSM. We first note that the 
energy-momentum tensor in the CQSM formally takes the same form as that of QCD, since 
the effective pion degrees of freedom contained in it is not an independent fields of quarks. 
The leading order contribution to the isoscalar d\ comes from the zeroth order term in the 
collective rotational velocity of the soliton |16j.[17j. so that it is given in the following form : 

d u+d = _ Nc (n\r 2 [a ■ r'p' ■ r - ~a-*p]\n). (18) 
2 3 

ra<0 



The symbol J2n<o denotes the summation over the occupied states (i.e., the discrete valence 
level (n = 0) plus the negative-energy Dirac-sea orbitals (n < 0)) in the hedgehog mean field. 
In the above equation, = ~ (~p + p") with "p and j? being the momentum operators 
respectively acting on the initial and final state wave functions. After some manipulation, 
di +d can be transformed into a form, which is convenient for numerical calculation : 

d u +d = _^_ MnNc ( n \r 2 \[Y 2 {r) x V]« x a] (0) \ n). (19) 

It is an easy exercise that this precisely coincides with the following expression 

d u + d = — -N c M N (n | 7° 7 3 {V 3 , |a;| 2 P 2 (cos^) } |n), (20) 

n<0 

derived by Schweitzer et al. in the same model by starting with the expression for the nucleon 
unpolarized GPD H u+d (x,!;,t). First pointing out that the valence quark contribution to 
di +d vanishes identically, i.e. c?" +d (val) = 0, they estimated the contribution of the polarized 
Dirac sea by means of a kind of derivative-expansion type approximation. They thus find that 
d^ +d = d^ +d (sea.) ~ —9.46 at the model energy scale around 600 MeV. 

However, we find no reason why d\ +d (val) vanishes. In fact, the relevant operator appearing 
in Eq. (fT9l) is a positive parity operator with the total grand spin K being zero (here K = 
J + 1 t), and there is no selection rule which enforces its matrix element between the valence 
quark state with the quantum numbers K p = + to vanish. (Although not so obvious, the 
operator appearing in Eq. (1201) also contains the K = component.) Here, we shall calculate 
this discrete valence level contribution explicitly. We also try to evaluate the contribution of the 
deformed Dirac sea, without recourse to the derivative expansion type approximation. This is 
possible with use of the discretized momentum basis of Kahana,Ripka and Soni [18], [19]. In the 
present analysis, we use the self-consistent soliton solutions obtained in [21] within the double- 
subtraction Pauli-Villars regularization scheme [20]. The model in the chiral limit contains 
only one parameter M, i.e. the dynamical quark mass. Here, we use the value M = 400 MeV. 
Since the energy-momentum-tensor distribution carried by the quark fields is expected to be 
sensitive to the value of the pion mass, we shall also investigate the pion mass dependence of 
d™ +d . The effective model lagrangian, which incorporates the finite pion mass effects, is given 
in [20]. Self-consistent soliton solutions are prepared in [21] for several values of pion mass. As 
pointed out in that paper, favorable physical predictions of the model are obtained by using the 
value of M = 400 MeV and m„ = 100 MeV, since this set gives a self-consistent solution close 
to the phenomenologically successful one obtained with M = 375 MeV and m n = MeV in the 
single-subtraction Pauli-Villars regularization scheme in the previous studies of nucleon parton 
distribution functions [22] -[28]. We first show the theoretical prediction for d\ +d corresponding 
to this favorable parameter set. It gives 

d u +d = d u+d( var) + rf «+d( sea ) ; (21) 
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with 



^ +d (val) ~ 0.66, < +d (sea) 



5.51. 



(22) 



We confirm that the dominant contribution to comes from the quarks in the negative- 
energy Dirac-sea orbitals. This deformed Dirac-sea contribution is large and negative. How- 
ever, we also find that the quarks in the discrete valence level gives nonzero and positive con- 
tribution, which partially cancels the Dirac-sea contribution. As pointed out in [5]? the large 
and negative prediction for the D-term coefficient is a characteristic feature of the CQSM, 
which maximally incorporates the spontaneous breaking of the chiral symmetry. In fact, the 
dominant contribution from the polarized Dirac sea is also viewed as simulating the t-channel 
exchange of two pions with the quantum numbers J PG = ++ ,2 ++ , • • • [6]. The final model 
predictions for d± +d obtained as the sum of the valence and the Dirac-sea contribution is 



which is a little smaller than the prediction d% +d ~ —8.0 obtained from the numerically evalu- 
ated H u+d (x,£,t) [TTj.flO] and the prediction d^ +d ~ —9.46 obtained based on the derivative- 
expansion type approximation with neglect of the valence level contribution in the same model 
[j~2] . In view of the difference of the soliton profile functions used in all these analyses (note 
that our result corresponds to d^ +d ~ —6.2 in the chiral limit), the qualitative agreement 
is encouraging, and it confirms the unique feature of the CQSM, which takes account not 
only of three valence quarks but also infinitely many Dirac-sea quarks in the mean potential. 
This is clear, if one compares the above predictions of the CQSM with that of the MIT bag 
model, d^ +d (MTT) ~ —0.716. We also recall that the lattice QCD simulation performed by the 
QCDSF collaboration in the heavy pion region around ~ 800 MeV [13] gives fairly small 
number : < +d (QCDSF) = § C 2 u+d (0) ~ -0.25 ± 0.13. Although the lattice QCD prediction 
quoted here corresponds to the energy scale around Q 2 ~ 4(GeV) 2 , there seems to be more 
difference in magnitude than explained by its scale dependence. 

The pion cloud interpretation of the Dirac-sea contribution to d^ +d in the CQSM may 
also be confirmed by investigating its pion mass dependence. We show in Fig.l the CQSM 
prediction for d\ +d in dependence of the pion mass. One sees that the magnitude of d^ +d 
rapidly decreases as increases, showing a tendency to match the very small lattice QCD 
prediction obtained in the heavy pion region at least qualitatively. Since the valence quark 
contribution to d™ +d is less sensitive to the variation of the pion mass, it can be interpreted as 
the reduction of the pion cloud effects as the model parameter m n increases. 

Now, we turn to the discussion of the isovector d\- It was emphasized in [5] that the 
isovector D-term is suppressed relative to the isoscalar one by a factor of l/N c , so that it is 
negligible in the large N c limit. Since N c = 3 in reality, however, it is not self-evident whether 
it is in fact numerically small or not. Here, we try to evaluate this subleading term in the 




-4.85, 



(23) 
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Figure 1: The CQSM prediction for d± +d in dependence of the pion mass. 

1/N C expansion in an explicit manner. In conformity with the observation above, d^~ d in the 
CQSM survives only at the 1st order in the collective rotational velocity of the soliton, so that 
the answer is given in the following form, i.e. as a double sum over the single quark orbitals 
in the hedgehog mean field : 

v u m>0,n<0 m n 

x (m\\r 2 [[y 2 (f) x Vp x a] (0) r\\n). (24) 

Let us first show the prediction obtained with the favorable set of the model parameters, 
i.e. M = 400 MeV and m n = 100 MeV. This gives 

d u ~ d = <- d (val) + <- d (sea) ~ 0.33 + 0.01 ~ 0.34. (25) 

We find that, in contrast to the isoscalar case, the contribution of the discrete valence level 
is dominant, while that of the deformed Dirac-sea is almost negligible. We also confirm that 
net di~ d is much smaller in magnitude than d^ +d , in conformity with the expectation based 
on the large N c counting. Still, our explicit calculation shows that the difference of d\ and df 
takes nonzero positive value owing to the presence of the discrete valence level contribution. 
Also interesting is how the isovector d± depends on the pion mass. We show in Fig. 2 the 
CQSM prediction for di~ d in dependence of the pion mass. Contrary to the isoscalar case, the 
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Figure 2: The CQSM prediction for d\~ d in dependence of the pion mass. 

magnitude of is an increasing function of m^. This peculiar behavior of d^~ d resembles 
the pion mass dependence of the generalized form factor B2Q d (t) at t = 0, or equivalently the 
isovector gravito anomalous magnetic moment of the nucleon, so that it may have a similar 
origin. (See Fig. 5 of [21], and the explanation around there.) 

In summary, we have carried out a theoretical analysis of the most important constants 
that characterize the nucleon D-term, i.e. the first coefficients di of its Gegenbauer expan- 
sion, within the framework of the CQSM, without recourse to the derivative-expansion type 
approximation used in the previous studies. We gave predictions not only for the leading 
isoscalar part of di but also for the subleading isovector part of d\ in the 1/N C expansion. 
Our treatment makes it possible to estimate the contribution of the discrete valence level and 
that of the negative-energy Dirac-sea levels separately. We found that, as for the isoscalar 
, the contribution of the deformed Dirac-sea is large and negative and dominates over 
the small positive contribution from the discrete valence level. This reconfirms the unique 
feature of the CQSM, which takes good account of the effects of the pion cloud generated 
by the spontaneous chiral symmetry breaking of the QCD vacuum. The predicted value of 
^u+d ^ _(4.9 6.2) at the model energy scale around 600 MeV qualitatively supports the 
previous estimate d^ +d ~ — (8.0 ~ 9.6) obtained in the same model based on the derivative- 
expansion type approximation. We have also found that d%~ ~ 0.34, which we confirm is 
much smaller than d" + but cannot be completely neglected. We hope that the theoretical 
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analysis carried out here will give useful constraints on the future analyses of high-energy 
DVCS and DVMP processes using the double distribution parameterization of the nucleon 
GPDs supplemented with the D-term. 
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Note added 

After submission of the present paper, two papers [29] and [30] appeared where the constant 
d\ +d (in addition to some other nucleon form factors of the energy momentum tensor) was 
calculated in a different manner within the same model with different regularization scheme. 
In the footnote 4 of |29j, it was pointed out that the statement di +d (val) = made in [12] is 
incorrect, which is consistent with our observation in the present paper. 
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